Astronomy & Astrophysics manuscript no. a773astroph 


© ESO 2008 


February 5, 2008 





Internal dynamics of the radio halo cluster Abell 773: 
a multiwavelength analysis 

R. Barrena 1 , W. Boschin 2 - 3 , M. Girardi 3 ' 4 , and M. Spolaor 3 ' 5 

1 Instituto de Astrofisica de Canarias, C/Via Lactea s/n, E-38205 La Laguna, Tenerife, Canary Islands, Spain 

2 Fundacion Galileo Galilei - INAF, C/Alvarez de Abreu 70, E-38700 Santa Cruz de La Palma, Canary Islands, Spain 

3 Dipartimento di Astronomia, Universita degli Studi di Trieste, via Tiepolo 11, 1-34131 Trieste, Italy 

4 INAF - Osservatorio Astronomico di Trieste, via Tiepolo 11, 1-3413 1 Trieste, Italy 

5 Centre for Astrophysics & Supercomputing, Swinburne University, Hawthorn, VIC 3122, Australia 

Received / accepted 

ABSTRACT 

Aims. To conduct an intensive study of the rich, X-ray luminous galaxy cluster Abell 773 at z = 0.22 containing a diffuse radio halo 
to determine its dynamical status. 

Methods. Our analysis is based on new spectroscopic data obtained at the TNG telescope for 107 galaxies, 37 spectra recovered from 
the CFHT archive and new photometric data obtained at the Isaac Newton Telescope. We use statistical tools to select 100 cluster 
members (out to ~1.8 fcl Mpc from the cluster centre), to analyse the kinematics of cluster galaxies and to determine the cluster 
structure. Our analysis is also performed by using X-ray data stored in the Chandra archive. 

Results. The 2D distribution shows two significant peaks separated by ~2' in the EW direction with the main western one closely 
located at the position of the two dominant galaxies and the X-ray peak. The velocity distribution of cluster galaxies shows two peaks 
at v ~65000 and ~67500 km s -1 , corresponding to the velocities of the two dominant galaxies. The low velocity structure has a high 
velocity dispersion - cr v = 800 - 1 100 km s _I - and its galaxies are centred on the western 2D peak. The high velocity structure has 
intermediate velocity dispersion - <x v ~ 500 km s _I - and is characterized by a complex 2D structure with a component centred on 
the western 2D peak and a component centred on the eastern 2D peak, these components probably being two independent groups. We 
estimate a cluster mass within 1 fcl Mpc of 6-1 1 xlO 14 M Q . Our analysis of Chandra data shows the presence of two very close 
peaks in the core and the elongation of the X-ray emission in the NEE-SWW direction. 

Conclusions. Our results suggest we are looking at probably two groups in an advanced stage of merging with a main cluster and 
having an impact velocity Av r f ~ 2300 km s . In particular, the radio halo seems to be related to the merger of the eastern group. 

Key words. Galaxies: clusters: general - Galaxies: clusters: individual: Abell 773 - Galaxies: distances and redshifts - Cosmology: 
observations 



1. Introduction 

Clusters of galaxies are recognized to be not simple relaxed 
structures, but rather as evolving via merging processes in a 
hierarchical fashion from poor groups to rich clusters. Much 
progress has been made in recent years in the observations of the 
signatures of merging processes (see Feretti et al. 2002 for a gen- 
eral review). The presence of substructure, which is indicative of 
a cluster at an early phase of the process of dynamical relaxation 
or of secondary infall of clumps into already virialized clusters 
occurs in about 50% of clusters as shown by optical and X-ray 
data ( e.g. Geller & Bee rs [T982l Mohr et al. [T996l G irardi et al. 
[19971 Kries sler & Beers [1997l J ones & Forman [1999l Schuecker 
et al. 120011 Burgett et al. |2004l . 

A very interesting aspect of these investigations is the possi- 
ble connection of cluster mergers with the presence of extended 
diffuse radio sources, halos and relics. They are rare, large (up 
to ~ 1-2 fcl Mpc), amorphous cluster sources of uncertain ori- 
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gin and generally steep radio spectra (Hanisch 1982; see also 
Giovannini & Feretti 2002] for a recent review). They appear to 
be associated with very rich clusters that have undergone recent 
mergers, and it has therefore been suggested by various authors 
that cluster halos/relics are related to recent merger activity (e.g. 
Tribble[[993l Burns et al. [19941 Feretti |T999>. 

The synchrotron radio emission of halos and relics demon- 
strates the existence of large scale cluster magnetic fields, of the 
order of 0. 1-1 juG, and of widespread relativistic particles of en- 
ergy density 10~ 14 - 10~ 13 erg cm . The difficulty in explaining 
radio halos arises from the combination of their large size and the 
short synchrotron lifetime of relativistic electrons. The expected 
diffusion velocity of the electron population is of the order of the 
Alfven speed (~100 km s~') making it difficult for the electrons 
to diffuse over a megaparsec-scale region within their radiative 
lifetime. Therefore, one needs a mechanism by which the rel- 
ativistic electron population can be transported over large dis- 
tances in a short time, or a mechanism by which the local elec- 
tron population is reaccelerated and the local magnetic fields are 
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amplified over an extended region. The cluster-cluster merger 
can potentially supply both mechanisms (e.g. Giovannini et al. 
119931 Bu rns et a ll 19941 Rottgering et al. 119941 see also Feretti et 
al. 2002 ; Sarazin 2002 for reviews). However, the question is still 
debated since the diffuse radio sources are quite uncommon and 
only recently we have been able to study these phenomena on the 
basis of sufficient statistics (a few dozen clusters up to z ~ 0.3, 
e.g. Giovannini et al. 119991 see also Giovannini & Feretti 2002 ; 
Feretti|2005j. 

Growing evidence for the connection between diffuse ra- 
dio emission and cluster merging is based on X-ray data (e.g. 
Bohringer & Schuecker 2002J Buote |2002| l. Studies based on 
a large number of clusters have found a significant relation be- 
tween the radio and the X-ray surface brightness (Govoni et al. 
2001a, 2001b) and between the presence of radio halos/relics 
and irregular and bimodal X-ray surface brightness distribu- 
tion (Schuecker et al. 2001). New unprecedent insights into 
merging processes in radio clusters are offered by Chandra and 
XMM-Newton o bservat ions (e.g. Ma rkevitc h & Vikhlinin l200U 
Markevitch et al. [20021 Fujita et al. [20041 Henry et al. [20041 
Kempner & David [2004b . 

Optical data are a powerful way to investigate the presence 
and the dynamics of cluster mergers, too (e.g. Girardi & Biviano 
2002 ). The spatial and kinematic analysis of member galaxies al- 
low us to detect and measure the amount of substructure, to iden- 
tify and analyse possible pre-merging clumps or merger rem- 
nants. This optical information is complementary to X-ray in- 
formation since galaxies and ICM react on different time scales 
during a merger (see, for example, numerical simulations by 
Roettiger et al. I1997I ). Unfortunately, to date optical data are 
lacking or are poorly exploited. The sparse literature contains 
a few individual clusters (e.g. Colless & Dunn 1996; Gomez et 
al. 120001 Barrena et al.[2002; Mercurio et al. [20031 Boschin et 
al. 120041 Boschin et al.[2006; Girardi et al. 2006). We have con- 
ducted an intensive study of Abell 773 (hereafter A773), which 
has a radio halo (Giovannini et al. 119991 see also Kempner & 
Sarazin |200T1 l. 

A773 is a rich, X-ray luminous and hot cluster at z ~ 0.22: 
Abell richness class = 2 (Abell et al. [19891), L x (0. 1-2.4 keV)= 
12.52 x 10 44 /i~ 2 erg s " 1 (Ebeling et al. [1996b , T x ~ 7-9 keV 
(Allen & Fabian TT99%l Rizzaet al. [19981 White[2000; Saunders 
et al. 120031 Govoni et al. 120051 Ota & Mitsuda[2005J. This clus- 
ter has a centre that is extremely optically rich and contains al- 
most only early-type galaxies (Dahle et al. 2002). A lot of evi- 
dence points towards the young dynamical status of A773. This 
cluster is elongated along the NE-SW direction and has a great 
amount of substructure in X-ray emission, as shown by ROSAT 
H RI data (Pierre & Starck [199%1 Rizza et al. [19981 Govoni et 
al. [2001b). There is a strong peak in the mass map, offset to the 
southwest from the optical centre by about 1' and the mass ap- 
pears to be elongated in a different way with respect to the light 
and number density distribution (Dahle et al. 120021) . The optical 
DSS image clearly shows two galaxy substructures, one in the 
centre of the X-ray emission and another at the eastern region 
(Govoni etal. l2004b . 

To date, few spectroscopic data have been reported in the lit- 
erature. Crawford et al. ( 1995) measured the redshift for the two 
brightest, equally dominant galaxies (z ~ 0.22). A few redshifts 
for radio galaxies in the cluster region are given by Morrison et 
al. (12003). We have recently carried out spectroscopic observa- 
tions with the TNG giving new redshift data for 107 galaxies in 
the field of A773, as well as photometric observations at the INT. 
We recover additional spectroscopic information from reduction 
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Fig. 1. 7?-band image of A773 (west at the top and north to the 
left) taken with the WFC camera of the INT. Galaxies with suc- 
cessful velocity measurements are labelled as in Table[3] Circles 
and boxes indicate cluster members and non-member galaxies, 
respectively. 

of CFHT archive data. Our present analysis is based on these 
optical data and X-ray Chandra archival data. 

This paper is organized as follows. We present the new op- 
tical data in Sect. 2 and the relevant analyses in Sect. 3. Our 
analysis of Chandra X-ray data is shown in Section 4. We dis- 
cuss the dynamical status of A773 in Sect. 5 and summarize our 
results in Section 6. 

Unless otherwise stated, we give errors at the 68% confi- 
dence level (hereafter c.l.). Throughout the paper, we assume a 
flat cosmology with Q m = 0.3, Oa = 0.7 and Ho - 70 hjo 
km s" 1 Mpc" 1 . For this cosmological model, 1' corresponds to 
213 kpc at the cluster redshift. 

2. Data sample 

2.1. Spectroscopy 

Multi-object spectroscopic observations of A773 were carried 
out at the TNG in December 2003 during the programme of 
proposal AOT8/CAT-G6. We used DOLORES/MOS with the 
LR-B Grism 1, giving a dispersion of 187 A/mm, and the Loral 
2048 x 2048 CCD (pixel size 15 //m). This combination of grat- 
ing and detector results in a dispersion of 2.8 A/pix. We have 
taken four MOS masks for a total of 137 slits. We acquired two 
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Fig. 2. /?-band image of the cluster A773 with, superimposed, 
the contour levels of the Chandra archival image ID 5006 (blue; 
photons in the energy range 0.3-7 keV) and the contour levels 
of a VLA radio image (red) at 1.4 GHz (Govoni et al. I2001bl >. 
Green ellipses identify structures detected by Wavdetect. The di- 
rection of A782 is also shown (see Sect. 5.4). North is at the top 
and east to the left. 

exposures of 1800 s for one mask and three exposures of 1800 
s for three masks. Wavelength calibration was performed using 
helium-argon lamps. 

Reduction of spectroscopic data was carried out with the 
IRAF package. Q 

Radial velocities were determined using the cross- 
correlation technique (Tonry & Davis 1979) implemented 
in the RVSAO package (developed at the Smithsonian 
Astrophysical Observatory Telescope Data Center). Each 
spectrum was correlated against six templates for a variety 
of galaxy spectral types: E, SO, Sa, Sb, Sc and Ir (Kennicutt 
119921 ). The template producing the highest value of % i.e. the 
parameter given by RVSAO and related to the signal-to-noise of 
the correlation peak, was chosen. Moreover, all the spectra and 
their best correlation functions were examined visually to verify 
the redshift determination. In two cases (IDs 66 and 129; see 
Table [3} we took the EMS AO redshift as a reliable estimate of 
the redshift. We obtained redshifts for 107 galaxies. 

For six galaxies we obtained two redshift determinations of 
similar quality. This allows us to obtain a more rigorous esti- 
mate for the redshift errors since the nominal errors as given by 
the cross-correlation are known to be smaller than the true er- 
rors (e.g. Malumuth et al. |1992t Bardelli et al l 19941 Ellingson 
& Yee 119941 Quintana et al. I2000I I. For the six galaxies hav- 
ing two redshift determinations we fit the first determination vs. 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



the second one by using a straight line and considering errors in 
both coordinates (e.g. Press et al. |1992| l. The fitted line agrees 
with the one-to-one relation, but, when using the nominal cross- 
correlation errors, the small value of x 1 probability indicates a 
poor fit, suggesting the errors are underestimated. Only when 
nominal errors are multiplied by a factor of ~1.5 can the ob- 
served scatter be justified. We therefore assume hereafter that 
true errors are larger than nominal cross-correlation errors by a 
factor 1.5. For the six galaxies we used the average of the two 
redshift determinations and the corresponding error. 

We added to our data observations stored in the CFHT 
archive (proposal ID: 01 AF37): 1 MOS mask of three 1800 s ex- 
posures with 52 slits. We reduced the spectra with the same pro- 
cedure adopted for the TNG observations. We succeeded in ob- 
taining spectra for 37 galaxies, two of which are in common with 
TNG galaxies. The comparison of the two estimates agree well 
within the errors (czcfht = 44170+199 vs. cz T ng = 44195 ±94 
km s _1 , czcfht = 64304 + 113 vs. cz T ng = 64260 + 46 km s" 1 ), 
allowing us to combine the data. For the two galaxies in com- 
mon we used the average of the two redshift determinations and 
the corresponding error. 

Our spectroscopic catalogue consists of 142 galaxies in a re- 
gion of 18'xl3' centred on the position of the two close domi- 
nant galaxies (IDs 59 and 60, hereafter Dl and D2; see Table[3]l. 
The median S/N and error on cz are 12 and 72 km s _1 , respec- 
tively. 

We use a conservative approach leading to a sparse spec- 
tral classification (75% of the sample, see Table [3). We follow 
the classification of Dressier et al. ( 119991 see also Poggianti et 
al. 119991 ). We define "e"-type galaxies as those showing active 
star formation indicated by the presence of [Oil], and, in par- 
ticular, "e(b)" galaxies as those showing an equivalent width 
EW([OII]) of < -40 A (probably starburst galaxies). Of galax- 
ies having S/N > 10 we define as "k"-type those having 

EW(H5)<3 A and no emission lines (probably passive galaxies); 
"k+a"- and "a+k"-type galaxies as those having 3 <EW(H5)<8 
A and EW(H<5)>8 A, respectively, and no emission lines (the 
"post-starbust" galaxies of Couch et al. 1994); "e(a)"-type galax- 
ies as those having emission lines and EW(H<5)>4 A ; "e(c)"- 
type galaxies as those having emission lines and EW(H5)<4 
A (probably spiral galaxies). Galaxy ID 89, which is the only 
one showing any emission lines ([OIII] in this case), but having 
the [Oil] line outside the observed spectral range, is classified 
as "e(a)". Out of 106 classified galaxies (94 with S /N > 10) 
we find 68 "passive" galaxies and 38 "active" galaxies (starbust 
and poststarbust; i.e. 6:1:8:9:5:9 for k+a:a+k:e(c):e(a):e(b):e, re- 
spectively). Of nine galaxies of type "e", i.e. those moderate 
emission line galaxies with low S/N, six are probably of type 
e(c) and three are probably of type e(a) (hereafter e(c?) and 
e(a?)). 

2.2. Photometry 

Our photometric observations were carried out with the Wide 
Field Camera (WFC), mounted at the prime focus of the 2.5 m 
INT (located at Roque de los Muchachos Observatory, La Palma, 
Spain). We observed A773 in December 2004 in photometric 
conditions with a seeing of about 1 .4" . 

The WFC consists of a four-CCD mosaic covering a 
33'x33' field of view, with only a 20% marginally vignetted 
area. We took nine exposures of 720 s in B H and 300 s in R H 
Harris filters (a total of 6480 s and 2700 s in each band) devel- 
oping a dithering pattern of nine positions. This observing mode 
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allowed us to build a "supersky" frame that was used to cor- 
rect our images for fringing patterns (Gullixson 1992). In addi- 
tion, the dithering helped us to clean cosmic rays and avoid gaps 
between the CCDs in the final images. The complete reduction 
process (including flat fielding, bias subtraction and bad-column 
elimination) yielded a final coadded image where the variation 
of the sky was lower than 0.4% in the whole frame. Another ef- 
fect associated with the wide field frames is the distortion of the 
field. In order to match the photometry of several filters, a good 
astrometric solution is needed to take into account these distor- 
tions. Using IRAF tasks and taking as a reference the USNO 
B1.0 catalogue, we were able to find an accurate astrometric so- 
lution (rms ~ 0.5") across the full frame. The photometric cal- 
ibration was performed using Landolt standard fields obtained 
during the observation. 

We finally identified galaxies in our Z?h and Ru images and 
measured their magnitudes with the SExtractor package (Bertin 
& Arnouts 1996) and AUTOMAG procedure. In a few cases (e.g. 
close companion galaxies, galaxies close to defects of the CCD) 
the standard SExtractor photometric procedure failed. In these 
cases we computed magnitudes by hand. This method consists 
in assuming a galaxy profile of a typical elliptical and scaling it 
to the maximum observed value. The integration of this profile 
gives us an estimate of the magnitude. This method is similar to 
PSF photometry, but assumes a galaxy profile, more appropriate 
in this case. 

We transformed all magnitudes into the Johnson- 
Cousins system (Johnson & Morgan 119531 Cousins 
[19761 We used B — Bu + 0.13 and R = R H 
as derived from the Harris filter characterization 
(http://www.ast.cam.ac.uk/~wfcsur/technical/photom/colours/) 
and assuming a B - V ~ 1.0 for E-type galaxies (Poggianti 
119971 ). As a final step, we estimated and corrected the galactic 
extinction As ~ 0.06, Ar ~ 0.04 from Burstein & Heiles's 
( 1982 ) reddening maps. 

We estimated that our photometric sample is complete down 
to R = 22 A (23.6) and B = 20 .4 (22.5) for S /N = 5 (3) within 
the observed field. 

We assigned magnitudes to 141 of the 142 galaxies of our 
spectroscopic catalogue. We measured redshifts for galaxies 
down to magnitude R ~ 21, but a high level of completeness 
is reached only for galaxies with magnitude R < 19.5 (~50% 
completeness). 

Table [3] (available in electronic format) lists the velocity 
catalogue (see also Fig. Q] and Fig. |2): identification number 
ID (Col. 1), identification label of each galaxy following IAU 
nomenclature rules (Col. 2); right ascension and declination, a 
and 6 (J2000, Col. 3); B magnitudes (Col. 4); R magnitudes 
(Col. 5); heliocentric radial velocities, v = czq (Col. 6) with 
errors, Av (Col. 7); spectral classification (Col. 8). 



3. Analysis of optical data 

3.1. Member selection and global properties 

To select cluster members from 142 galaxies having redshifts, 
we use the adaptive-kernel method (hereafter DEDICA, Pisani 
[19931 and [19961 se e also Fadda et al. [19961 Girardi et al. [19961 
Girardi & Mezzetti 2001). We find significant peaks in the veloc- 
ity distribution >99% c.l.. This procedure detects A773 as a two- 
peaked structure, populated by 100 galaxies with 62250 < v < 
68850 km s _1 , hereafter referred as cluster members (see Table[3] 
and Fig. [3j. Of the non-member galaxies, 24 and 18 are fore- 
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Fig. 3. Redshift galaxy distribution. The solid line histogram 
refers to galaxies assigned to the cluster according to the 
DEDICA reconstruction method. 



ground and background galaxies, respectively. In particular, 15 
foreground galaxies belong to a low density peak at z ~ 0.186. 

By applying the biweight estimator to cluster members 
(Beers et al. 1990), we compute a mean cluster redshift of 
<Z> = 0.2197+ 0.0005, i.e. (v) = 65876+140 km s _1 . We es- 
timate the LOS velocity dispersion, cr v , by using the biweight 
estimator and applying the cosmological correction and the stan- 
dard correction for velocity errors (Danese et al. 1980). We ob- 
tain <x v = 1394^g km s _1 , where errors are estimated through a 
bootstrap technique. 

Hereafter, we consider as cluster centre the position of 
the most luminous dominant galaxy [Dl, R.A.=09 h 17 m 53 ?26, 
Dec.=+51°43'36:'5 (J2000.0)], which is related to the richer 
subsystem as shown in the following sections. 



3.2. Bimodal velocity structure 

The DEDICA method assigns 63 galaxies with v < 66543 
km s _1 to a peak at cz = 65013 km s and 37 galaxies with 
v > 66653 km s -1 at a peak at cz = 67453 (z ~ 0.217 and 
z ~ 0.224, respectively, see Figure[4]i. According to the working 
definition of Girardi et al. (1996), the two peaks are not eas- 
ily separable since their overlap is very large, i.e. many galaxies 
(63/100) have a non-negligible probability of belonging to both 
peaks. 

The DEDICA method used above has the advantage of not 
requiring any a priori shape for the subsystem research. Kaye's 
Mixture Model (KMM, as implemented by Ashman et al. 1994 ), 
which fits a user-specified number of Gaussian distributions, 
gives a quite similar result. Of the cluster members we find that 
the KMM method fits a two-group partition by rejecting the sin- 
gle Gaussian at the 99.7% c.l., as obtained from the likelihood 
ratio test, assigning the 65 galaxies with v < 66659 km s _1 and 
the 35 galaxies v > 66845 km s _1 to two different groups. 
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Table 2. Results of the kinematical analysis 
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Fig. 4. Lower panel: rest-frame velocity vs. projected distance 
from the cluster centre - the location of the dominant galaxy 
Dl - of the 100 galaxies assigned to the cluster (Fig. [3}. Upper 
panel: velocity histogram of the 100 galaxies assigned to the 
cluster. Velocities of the two dominant galaxies are pointed out. 
The two-peaks detected by the DEDICA method and the two 
Gaussians corresponding to the KMM partition (see Tabled are 
also shown (faint and dashed lines, respectively). 



Finally, we note that the two dominant galaxies of the cluster 
are assigned to different groups: Dl lies very close to the low ve- 
locity peak and D2 lies very close to the peak of the high velocity 
peak (Fig.|4|. 

We present the main properties for low- and high-velocity 
structures as defined by both the DEDICA method (hereafter 
LV- and HV-groups, respectively) and the KMM method (here- 
after, KMM1 and KMM2, respectively). Table |2]lists the name 
of the sample (Col. 1); the number of assigned members, N g 
(Col. 2); the mean velocity and its error, < v > (Col. 3); the ve- 
locity dispersion and its bootstrap errors, cr v (Col. 4). The above 
properties are computed using the galaxies assigned to each in- 
dividual group. However, since there is a wide velocity range 
where galaxies have a non-zero probability of belonging to both 
groups, both DEDICA and KMM group membership assign- 
ment leads to an artificial truncation of the distributions being 
"deconvolved" (see the velocity thresholds above). This trun- 
cation may lead to an underestimate of velocity dispersion for 
the subclusters (see Bird 1994). Therefore, we also list the value 
of <x v as computed by KMM software using the galaxies of the 
whole sample weighted with their assignment probabilities to 
the respective sample after having applied the cosmological cor- 
rection and the standard correction for velocity errors (Danese et 
al. [19801 

Since the velocity difference between the two groups is prob- 
ably of a dynamical nature (see Sect. 5) we expect that the mag- 
nitude incompleteness of our sample affects both in a similar 
way, thus having a small influence on their detection and kine- 
matic properties. To verify this issue, of 142 galaxies having 
redshifts we select samples of 98 (62) galaxies brighter than 



Sample 


N 8 


< v > a 






km s 


Whole system 


100 


65876 ± 140 


LV - group 


63 


64875 ± 105 


HV - group 


37 


67664 + 86 


KMM1 


65 


64924 ± 107 


KMM2 


35 


67717 + 83 


W- sub 


14 


65573 ± 366 


E - sub 


5 


67469 ± 86 


LV - W - sub 


8 


64415 + 198 


HV - W - sub 


6 


67187 ±262 


LV - very red 


28 


64962 + 100 


HV - very red 


15 


67579 ± 138 



km s 1 



1394 68 

825+if 
515-1 

858!^(1122) c 

484^«(548) c 
j 299+220 

473 + ff 

493+Jf 

539* 

518! 

506+f 



For comparison v of Dl and D2 dominant galaxies are 65131 and 
67765 km s"'. 

We use the biweight and the gapper estimators by Beers et al. 
( 1990) for samples with N 5 > 15 and with N 5 < 15 galaxies, re- 
spectively (see also Girardi et al. 1993 1. 

Between parenthesis we list (standard) velocity dispersions as com- 
puted by using the galaxies of the whole system weighted with their 
assignment probabilities to the respective sample rather then only 
the galaxies assigned to the samples (see text). 



19.5 (19) mags for which the completeness is at least ~ 50%. 
In both cases the KMM method detects the two peaks at simi- 
lar mean velocities. The estimates of velocity dispersions agree 
within 1—1.5 sigma with the values obtained in the whole sam- 
ple (i.e. cr v = 750+g 6 , 671^1° km s" 1 and <r v = 514+j® 474+ 
km s _1 for the main and secondary peaks, respectively - to be 
compared with KMM1 and KMM2 values). 



3.3. 2D galaxy distribution 

By applying the DEDICA method to the 2D distribution of 
A773 galaxy members we find two peaks. The position of 
the highest peak (hereafter W-peak) is very close to the loca- 
tion of the two dominant galaxies and of the peak of X-ray 
emission (see Sect. 4), while the secondary peak (hereafter E- 
peak) lies ~ 2' east [R.A.=09 h 17 m 53 !22, Dec.=+51°43'46 '.'0 
and R.A.=09 h 18 m 04 ! 15, Dec.=+51°43'40 '.'9 (J2000.0), respec- 
tively, see Fig. [5). 

Our spectroscopic data do not cover the entire cluster field 
and suffer from magnitude incompleteness. To overcome these 
limitations we recover our photometric catalogue by selecting 
likely members on the basis of the colour-magnitude relation 
(hereafter CMR), which indicates the early-type galaxy locus. To 
determine the CMR we fix the slope according to Lopez-Cruz 
et al. ( 120041 see their fig. 3) and apply the two-sigma-clipping 
fitting procedure to the cluster members obtaining B — R — 
3.507-0.072x7? (see Fig . [6j . From our photometric catalogue we 
consider galaxies (objects with SExtractor stellar index < 0.9) 
lying within 0.25 mag of the CMR. To avoid contamination by 
field galaxies we do not show results for galaxies fainter than 
21 mag (in the /?-band). The contour map for 582 likely cluster 
members having R < 21 shows how the main structure is elon- 
gated in the EW direction with two main peaks corresponding 
to those determined above (see Fig. |7). To check the effect of 
the gaps between the CCD chips of the WFC we also consider 
the distribution of 563 galaxies having R < 2\ and lying more 
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Fig. 5. Spatial distribution on the sky of spectroscopically con- 
firmed cluster members and the relative isodensity contour map. 
The contour map shows two main peaks: the main, western peak 
(W-peak) and the secondary, eastern peak (E-peak). The clus- 
ter members are denoted by circles (LV-galaxies) and squares 
(HV-galaxies) and the two dominant galaxies are shown by solid 
symbols. The plot is centred on the cluster centre defined as the 
position of the dominant galaxy Dl (solid circle). 



than 0.75 mag from the colour-magnitude relation, i.e. proba- 
bly mainly formed by non-member galaxies. The comparison of 
members and non-members in Fig. [7] confirms the reality of the 
above two peaks in the cluster structure. 

When analysing the 2D galaxy distribution for magnitude 
intervals the bimodal structure in the central cluster region is 
confirmed by the 141 galaxies brighter than 19 mag (in the R- 
band), while two samples of fainter galaxies (177 with 19 < R < 
20 and 262 with 20 < R < 21) show only one significant peak 
closer to the W-peak than to the E-peak (see Figure[8]). 

3.4. 3D structure of A773 

The existence of correlations between positions and velocities of 
cluster galaxies is a footprint of real substructures. 

In the case of A773 the standard approaches to determine 
these correlations fail (e.g. Boschin et al. 120061 Girardi et 
al. 12006). In fact, we find no significant velocity gradient, no 
significant substructure according to the method by Dressier & 
Schectman (I1988K and the result of the 3D KMM method is 
fully driven by the velocity-variable making us obtain again the 
two-group partition of Section 3.2. 

As an alternative approach, we analyse the 2D distributions 
of the two groups determined using the velocity distribution. 
We compare the 2D distribution of the galaxies assigned to the 
LV-group (hereafter LV-galaxies) with that of the galaxies as- 
signed to the HV-group (hereafter HV-galaxies) finding no sig- 
nificant difference according to the 2D Kolmogorov-Smirnov 
tests (2DKS test, see Fasano & Franceschini 1987 as imple- 
mented by Press et al. 119921 ). 

To go deeper into the question we also use the 2D DEDICA 
method. We find that the distribution of LV-galaxies shows a 
main, very significant peak located close to the W-peak, while 
the main 2D peak of the distribution of HV-galaxies is closely 
located near the E-peak, thus indicating a relation between ve- 
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Fig. 6. B - R vs. R diagram for galaxies with available spec- 
troscopy is shown by circles and small crosses (cluster and field 
members, respectively). The solid line gives the best-fit colour- 
magnitude relation as determined from member galaxies; the 
dashed lines are drawn at +0.25 mag from the CMR. According 
to our working definition in Sect. 3.5, cluster members are di- 
vided into "very red", "red" and "blue" galaxies (above the 
CMR, between the CMR and the dotted line and below the dot- 
ted line, respectively). 



locity structure and 2D structure detected in the above two sec- 
tions (see Fig. [9]). 

Moreover, the distribution of HV-galaxies also shows a sec- 
ondary peak closely located near the W-peak. This secondary 
peak is only 93% significant and much less dense with respect 
to the main one (only 30%), but contains a similar number of 
galaxies (seven vs. eight). This secondary peak is also found 
when we consider only the 22 galaxies with a velocity higher 
than the peak velocity of 67453 km s and therefore we exclude 
the possibility that the 2D secondary peak is induced by the dif- 
ficulties in member assignment between the two velocity peaks 
(see Sect. 3.2). The likely conclusion is that some HV-galaxies 
are really connected with the high velocity D2 galaxy located in 
the region of the W-peak. 

As a final approach, we analyse the velocity distributions of 
the galaxy populations located in the regions of the E- and W- 
peaks. Following Girardi et al. (120051 1 we analyse the profiles of 
mean velocity and velocity dispersion of galaxy systems corre- 
sponding to the E- and W-peaks to attempt an independent anal- 
ysis. 

In the case of the E-peak the integral cr v profile shows a sharp 
increase starting from the peak centre (see Fig.lTOl. This is prob- 
ably due to the contamination of a galaxy clump with a different 
mean velocity (see also Girardi et al. 1996 in the case of A3391- 
A3395). To avoid possible problems of contamination we decide 
to consider a sample E-sub formed by only the first five galaxies 
closest to the E-peak. In the case of the W-peak the value of <x v 
computed using the first five galaxies is enough high - ~ 1200 
km s _1 - and the integral <x v profile shows an almost stable be- 
haviour out to ~0.2 h~Q Mpc , where it starts to increase (see 
Fig.fTTTi. We decided to consider a sample W-sub formed by the 
first fourteen galaxies closest to the W-peak. 
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Fig. 7. Spatial distribution on the sky and relative isodensity con- 
tour map of 582 likely cluster members (according to the colour- 
magnitude relation) with R < 21, obtained with the DEDICA 
method (black lines). For comparison we show the contour map 
of the 563 likely non-cluster members (grey lines). The plot is 
centred on the cluster centre. The circle indicates the 5' central 
region. 
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Fig. 8. Isodensity contour map (zoomed within 1 }cl Mpc from 
the cluster centre) of likely cluster members of different mag- 
nitude intervals. For the sake of clarity we only show those for 
R < 19 and 19 < R < 20, indicated by black solid lines and 
dashed lines. We also indicate the position of the X-ray cen- 
tre (large triangle), and the E- and W-peaks as recovered from 
galaxies having redshifts (small dots, see also Fig. |5j. 



Table [2] lists kinematic properties for E- and W-subs and 
shows that these groups have similar properties of the HV- 
subcluster and of the whole cluster, respectively. When compar- 
ing the kinematics of E-sub and W-sub we find that their velocity 
distributions are different at the 97.8% c.l. according to the KS 
test (e.g. Ledermann 1982, as implemented by Press et al. 1 19921 
for small samples), and that the mean velocity and the velocity 
dispersion differ at the 99.6% c.l and 93.5% c.l. according to 
the means test and F test, respectively (e.g. Press et al. 11 9921 ) . 
Moreover, W-sub has a complex structure since its velocity dis- 
tribution shows two peaks according to DEDICA (see the small 
upper-right panel of Fig. [TT). The low velocity peak (LV- W-sub) 
contains the Dl galaxy and resembles the mean velocity of LV- 
subcluster. The high velocity peak (HV- W-sub) contains the D2 
galaxy and resembles the mean velocity of HV-subcluster. 

3.5. Further insights using galaxy populations 

In the case of A773, strongly substructured, the analysis of red 
and/or bright galaxies might be useful to trace the cores of 
the subclusters (e.g. Lubin et al. 120001 Boschin et al. 120041 . 
Moreover, active star forming galaxies might be related to the 
cluster-cluster merging phenomena (e.g. Bekki 1999). 

We divide the sample (99 galaxies having magnitudes) into 
low- and high-luminosity subsamples, each with 49 galaxies, by 
using the median R magnitude = 19.08. The two subsamples do 
not differ in their velocity distribution according to the KS test 
and both show a bimodal velocity distribution according to the 
DEDICA method. In agreement with the photometric sample re- 
sults, luminous galaxies clearly show the presence of both E and 
W peaks (see Figure [T2l. 

Butcher & Oemler (1984) define as "blue" those galaxies 
with rest-frame B - V colours at least 0.2 mag bluer than that 
of the CMR, this translates to 0.447 mag for B - R colours at 




X (arcmin) 

Fig. 9. Comparison of isodensity contour maps recovered from 
the galaxies assigned to the low velocity peak (black solid line) 
and from galaxies assigned to the high velocity peak (light grey 
line). The position of peaks W and E are indicated by solid cir- 
cles. 



Z ~ 0.21 (see Mercurio et al. 120041 and Haines et al. [20041 for 
A209), and a comparable result is obtained by using the ta- 
bles of Fukugita et al. (1995). Therefore, we hereafter define 
(B - R) C0lT = (B -R)- (3.507 - 0.072 x R) and "blue" those 
galaxies with (B - R) con < -0.45, i.e. those with observed B — V 
colours at least 0.45 mag bluer than that of the CMR. To analyse 
the colour segregation we also divide other galaxies into "red" 
[-0.45 < (B-R) con < 0] and "very red" [(B-R) con - > 0]. We ob- 
tain three subsamples of 43 very red, 39 red, and 17 blue galax- 
ies (see Fig. |6). The subsample of the very red galaxies shows 
two distinct peaks in the velocity distribution (LV-very red and 
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Fig. 10. Kinematic profiles of the cluster obtained assuming the 
E-peak as centre. The vertical line indicates the region which 
encloses the galaxies of E-sub (see text). Differential (big cir- 
cles) and integral (small points) mean velocity and LOS velocity 
dispersion profiles are shown in top and bottom panels, respec- 
tively. For the differential profiles are plotted the values obtained 
binning galaxies 10 by 10. For the integral profiles, the mean 
and dispersion at a given (projected) radius from the position 
of the E-peak is estimated by considering all galaxies within 
that radius: the first point refers to the first five galaxies and the 
last point refer to all galaxies contained in the cluster. The error 
bands at the 68% c.l. are also shown. 



HV-very red, see Table|2]i. When analysing the 2D galaxy distri- 
bution of each subsample through the DEDICA method we find 
that the very red galaxies show a significant peak at the position 
of the W-peak - this is true for both LV-very red and HV-very 
red subsamples separately - while red galaxies show a signifi- 
cant peak at the position of the E-peak (see Fig. [T3l ). The blue 
galaxies are instead preferentially located towards western re- 
gions of the cluster. 

Finally, we consider galaxies of different spectral types. Of 
76 classified member galaxies, we find 60 "passive" k-galaxies 
and 16 "active" galaxies, of which 10 are emission line galaxies 
(i.e. 5:1:1:3:1:5 for k+a:a+k:e(c)/e(c?):e(a)/e(a?):e(b):e, respec- 
tively). Both passive and emission line galaxy populations show 
the characteristic bimodal velocity distribution and have simi- 
lar mean velocity and velocity dispersion, quite comparable to 
that of the whole sample. As for the 2D distribution, the passive 
population shows the characteristic bimodal distribution with the 
E- and W-peaks, while the active populations avoid these re- 
gions of high density. In particular, Fig. [T4l shows that emission 
line galaxies are preferentially located towards south-western re- 
gions of the cluster. Indeed, we find that passive galaxies differ 
from the emission-line galaxies at the 98.0% c.l., respectively, 
according to the 2DKS test and despite some categories drawn 
in Fig. [14] include small numbers of objects and make this com- 
parison difficult. 



Fig. 11. The same as in Fig. [10] but referring to the W-peak. The 
small upper-right panel shows the velocity galaxy density (in 
normalized units), as provided by DEDICA for the galaxies of 
W-sub, i.e. those galaxies contained within the radius indicated 
by the vertical line. 
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Fig. 12. Spatial distribution on the sky of 100 cluster members. 
More and less luminous galaxies are denoted by circles and 
crosses, respectively. The contour map recovered from the most 
luminous galaxies is also shown. The position of peaks W and E 
are indicated by solid, small circles. 



4. X-ray analysis 

X-ray analysis of A773 has already been performed in the past, 
both based o n ROSA T (Pierre & Starck fTWgl Rizza et al.[l998 
Govoni et al. |2001bl > and Chandra data (Govoni et al. 120041 ). In 
particular, Govoni et al. (2004) analysed two Chandra ACIS-I 
observations (exposures ID #533 and #3588). The two pointings 
have a clean total exposure time of 18.7 ks (see their table 1). 
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Fig. 13. Spatial distribution on the sky of 100 cluster members. 
Very red, red and and blue galaxies are denoted by crosses, cir- 
cles and solid squares, respectively. The contour map recovered 
from the very red galaxies is also shown. The position of peaks 
W and E are indicated by solid small circles. 
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Fig. 14. Spatial distribution on the sky of 100 cluster members. 
Large symbols denote spectroscopically classified galaxies: cir- 
cles, squares, solid triangles, solid squares and solid pentagons 
indicate k,k+a/a+k,e(c)/e(c?),e(a)/e(a?) and e(b), respectively. 
Small dots indicate unclassified members. 



More recently, Chandra has again observed A773 (exposure ID 
#5006). This pointing has an exposure time of 20.1 ks. In this 
section we propose an analysis of this Chandra ACIS-I observa- 
tion. 

Data reduction is performed using the package CIAC0 
(Chandra Interactive Analysis of Observations) on chip 13 (field 
of view ~ 8.5' x 8.5')- First, we remove events from the level 2 
event list with a status not equal to zero and with grades one, five 
and seven. We then select all events with energy between 0.3 and 
10 keV. In addition, we clean bad offsets and examine the data, 
filtering out bad columns and removing times when the count 

2 CIAO is freely available at http://asc.harvard.edu/ciao/ 



rate exceeds three standard deviations from the mean count rate 
per 3.3 s interval. We then clean the 13 chip for flickering pixels, 
i.e. times where a pixel has events in two sequential 3.3 s in- 
tervals. The resulting exposure time for the reduced data is 19.6 
ks. 

In Fig. |2] we plot an /?-band image of the cluster with su- 
perimposed the X-ray contour levels of the Chandra image. The 
shape of the cluster appears to be moderately elliptical. By using 
the CIAO package Sherpa we fit an elliptical 2D Beta model to 
the X-ray photon distribution to quantify the departure from the 
spherical shape. The model is defined as follows: 



/(*,y)=/(r)=A/[l+(r/ro)T, 



(1) 



where the radial coordinate r is defined as r(x,y) - [X 2 (l - 
e) 2 + F 2 ] 1/2 /(l - e), X = (x - x ) cos 6 + (y - y ) sin 9 and 
Y = (y-yo) cos 9-{x-xq) sin 6. Here x and y are physical pixel 
coordinates on chip 13. The best fit centroid position is located 
very close to the dominant galaxy Dl, i.e. R.A.=09 h 17 m 52 ?87, 
Dec.=+51°43'39:'2 (J2000.0). The best-fit core radius, the el- 
lipticity and the position angle are ro = 47+3 arcsec (167+11 
h'l kpc), e = 0.23+0.02 and PA=79.7±2.8 degrees (measured 
counterclockwise from the north), respectively. 

Isophotes in Fig. [2] suggest that the core of A773 is substruc- 
tured with two distinct maxima in the surface brightness distribu- 
tion. The presence of these two peaks is confirmed by a wavelet 
analysis of the image performed by the task CIAO/Wavdetect. 
The task was run on several scales to search for substructures 
with different sizes. The significance thresholcQ was set at 10~ 6 . 
The results are shown in Fig. [2] Green ellipses represent two sig- 
nificant surface brightness peaks found by Wavdetect in the core 
of the cluster. They are located at R.A. 09 h 17 m 53.4 s and Dec. 
+51°43'44", and at R.A. 09 h 17 m 51.4 s and Dec. +51°43'32", 
respectively. This finding confirms the results of Pierre & Stark 
( 1998) and Rizza et al. (1998), who first noted the complex struc- 
ture of the core from a 16 ks ROSAT/HRI exposure. 

For the spectral analysis of the cluster X-ray photons, we 
compute a global estimate of the ICM temperature. The temper- 
ature is computed from the spectrum of the cluster within a cir- 
cular aperture of 2' radius around the cluster centre. Fixing the 
absorbing galactic hydrogen column density at 1.44xl0 20 cm" 2 , 
computed from the HI maps by Dickey & Lockman ( 1990), we 
fit a Raymond-Smith ( 1 19771 1 spectrum using the CIAO package 
Sherpa with a^ 2 statistics. Considering the energy range 0.8-10 
keV, we find a best fit temperature of Tx = 7.8 + 0.4 keV and 
a metal abundance of 0.34+0.07 solar, in agreement with the 
results of Govoni et al. (120041 ) . 



5. Discussion 

The global LOS velocity dispersion cr v = 1394^ km s _1 is 
considerably higher than the average X-ray temperature 
Tx =7.8±0.4 keV from our analysis of Chandra data when 
assuming the equipartition of energy density between gas and 
galaxies (i.e. y6 sp ec = 1.51+qJq to be compared with y6 spec = 10). 
This suggests that the cluster is far from dynamical equilibrium 
and indeed we do find evidence that A773 shows two peaks both 
in the velocity space (Av ~ 2800 km s _1 ) and in the 2D projected 
position space (~ 2' between the E-peak and the W-peak). 



see § 1 1 . 1 of the CIAO Detect Manual (software release version 3.3, 
available at the WWW site http://cxc.harvard.edu/ciao/manuals.html I 

4 Apec = o~l/(kT/fim p ) with = 0.58 the mean molecular weight and 
m p the proton mass. 
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The main subsystem is detected as the LV-peak in the ve- 
locity distribution, has a velocity dispersion of <x v c =800-l 100 
km s - comparable to massive clusters - and its centre prob- 
ably coincides with the Dl galaxy, located close to the main 2D 
peak (W-peak) and the peak of X-ray emission. The distribution 
of very red galaxies is also peaked in the same place. Both the 
low velocity galaxies located close to the W-peak (LV-W-sub) 
and the low velocity, very red galaxies (LV-very red) show a <x v 
value of about 500 km s _1 , lower than the global sigma as often 
observed in the case of a core cluster, in particular when exam- 
ining passive galaxies (Adami et al. [1998 , Biviano & Katgert 
120041 . 

A secondary subsystem is detected as the HV-peak in the 
velocity distribution and has a moderate velocity dispersion of 
cr v g ~ 500 km s . Its spatial structure is rather complex, show- 
ing two galaxy concentrations, one at the location of the E-peak 
and the other at the location of the W-peak, i.e. close to the high 
velocity, D2 dominant galaxy (see Fig. [9] and the upper-right 
panel of Fig. fTTb . The absence of a concentration of very red 
galaxies and of any close X-ray peak, as well as the presence 
of only a minor DM peak (Dahle et al. 120021) . suggest that the 
E-peak is dynamically of minor importance. Analysing galaxies 
around the two concentrations (E-sub and HV-W-sub) we find 
similar kinematic properties; thus, with present data we cannot 
distinguish whether they are two independent galaxy groups or 
two components of the same group. In favour of the first hy- 
pothesis, we notice that both theoretical predictions and obser- 
vational data suggest that clusters are formed by numerous small 
infalling systems (see the case of A85-A87 complex in Durret 
et al. 119981 1. In this case, the similarity of the mean velocities 
of E-sub and HV-W-sub might reflect that their velocities are of 
dynamical origin, both groups strongly interacting with the same 
cluster potential and thus having a similar merging velocity. In 
the case of the second hypothesis, E-sub could be either a dense 
substructure of the HV-system or a part of the core of the HV- 
system dynamically decoupled from the dominant galaxy during 
the cluster merger. We notice that Adami et al. ( 2005 ) found a 
similar feature in the Coma cluster, where the N03-SW group 
has the same mean velocity as NGC 4889, one of the two domi- 
nant galaxies, even if it is spatially distant in 2D. 

5.1. A773 mass estimate 

On the basis of the above discussion we estimate the mass of 
A773, assuming that the system is formed by a main cluster with 
a velocity dispersion of <x v c =800-l 100 km s _1 and one or two 
groups with cr v g ~ 500 km s . Here we report our results in 
the case of one group, but notice that in the case of two groups, 
they have the same virial masses and the results of this section 
are still valid. 

Although the system is probably in a phase of strong inter- 
action (see also the following section) so that the subsystems are 
still detectable in the velocity distribution, let us assume having 
a good estimate of cr v of the subsystems before the merger. To 
apply the virial method we must assume that each subsystem is 
in dynamical equilibrium. This might be not true in the case of a 
possibly substructured group (see above). On the other hand, the 
mass estimate of the whole complex depends above all on the 
main subsystem, whose structure is much more regular. 

Following the prescriptions of Girardi & Mezzetti ( 120011 1, 
we assume for the radius of the quasi-virialized region Ryir.c = 
O.nxcrJH(z) = 2.3-3.2/17,, 1 Mpc andR vil , g = 0.17xcr v ///( z ) ~ 
1.4 h~Q Mpc for the cluster and the group, respectively - see 



their eq. 1 after introducing the scaling with H(z) (see also eq. 8 
of Carlberg et al. 1 19971 for R2oo)- Therefore, our spectroscopic 
catalogue samples the whole group virialized region and more 
than half of the cluster one. 

One can compute the mass using the virial theorem (Limber 
& Mathews 1 19601 see also, for example, Girardi et al. 1998]! un- 
der the assumption that mass follows galaxy distribution: M = 
M sv ; r - SPT, where M sv ; r = 37r/2 • <XyR PV /G is the standard virial 
mass, Rpy a projected radius (equal two times the harmonic ra- 
dius), and SPT is the surface pressure term correction (The & 
White [T9861 . The estimate of cr v is generally robust when com- 
puted within a large cluster region (see Fig. 5 of Girardi et al. 
2006 and Fadda et al. |1996| for other examples), hence for each of 
the two subsystems we consider the corresponding above-global 
values cr v c and <x v g . The value of Rpy depends on the size of the 
region considered so that the computed mass increases (but not 
linearly) with the increasing region considered. Since the two 
subsystems are partially aligned along the LOS we avoid com- 
puting Rpy by using data for observed galaxies and use an al- 
ternative estimate that was shown to be good when Rpy is com- 
puted within R v ; r (see eq. 13 of Girardi et al. 1998). This al- 
ternative estimate is based on our knowledge of the galaxy dis- 
tribution and, in particular, a galaxy King-like distribution with 
parameters typical of nearby/medium-redshift clusters: a core 
radius R core = 1 /20 x R v i r and a slope-parameter {1^ = 0.8, 
i.e. the volume galaxy density at large radii as r~ 3 ^ r " = r~ 24 
(Girardi & Mezzetti 2001). For the whole virialized region we 
obtain R PV , C = 1.7-2.4 Mpc and R PV , g ~ 1.4 fc^J Mpc , 
where a 25% error is expected because typical, rather than in- 
dividual, galaxy distribution parameters are assumed. As for 
the SPT correction, we assume a 20% one, computed combin- 
ing data on many clusters (e.g. Carlberg et al. 119971 Girardi 
et al. 119981 1. This leads to virial masses for the subsystems of 
M c (< R vir , c = 2.3 - 3.2 /z^Mpc) = 1.0-2.5 xlO 15 h~\ M of 
M g (< R v ir,g = 1-4 /^Mpc) ~ 2.4 xlO 14 M Q with a typical 
error of ~30%. A reliable mass estimate of the whole system is 
then M = 1.2-2.7 xlO 15 hz^ M Q , in agreement with rich clus- 
ters reported in the literature (e.g. Girardi et al. 1998; Girardi & 
Mezzetti 120011 . 

To compare our result with the estimate obtained via grav- 
itational lensing we obtain a projected mass assuming that the 
cluster is described by a King-like mass distribution (see above), 
or a NFW profile where the mass-dependent concentration pa- 
rameter is taken from Navarro et al. (119971 1 and rescaled by the 
factor 1 + z (Bullock et al. [200T1 Dolag et al. 120041 . Moreover, 
we assume that the two subsystems have coincident centres, and 
that the subsystems themselves are extended only to one R v ; r . We 
obtain M proj (< R = 90")=(4.4-7 9)xl0 14 M Q in agreement 
with that found by Dahle et al. (2002) considering the quoted 
errors (Moahie ~ 3xl0 14 hzl M e with an upper error of the 50% 
after rescaling to our cosmological model). 

Using the same mass distributions we compute M(< R = 
1 Mpc)=(5.9-ll.l)xl0 1 4 M and, using r-band lu- 
minosity by Popesso et al. (2004), estimate M/L=( 150-280) 
hjo M Q /L Q . The upper value is typical of clusters of similar red- 
shift (Carlberg et al. [19961 . 

5.2. Merging scenario 

The X-ray emission shows two very close peaks along the NE- 
SW direction separated by ~0.5'. The main X-ray peak, the NE 
one, is closely located near - but is not coincident with - the Dl 
galaxy. The secondary, SW, one does not coincide with any lu- 
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minous galaxy or obvious galaxy substructure (see Fig.|2|. This 
lack of correlation between collisional and non-collisional com- 
ponents strongly suggests that the cluster is undergoing an ad- 
vanced phase of merging (e.g. Roettiger et al. |1997l l. Moreover, 
the direction of elongation of the isophotes of the X-ray emis- 
sion (NEE-SWW direction, PA~80 degrees) is similar to that 
indicated by the elongation of galaxy distribution (see, for ex- 
ample, contours for very red galaxies in Fig. [TBI , but quite dif- 
ferent from N-S one indicated by the two dominant galaxies and 
the elongation of the main DM clump (Dahle et al. 2002 ). It is 
worth noting that the literature reports several cases of clusters 
having two dominant galaxies with significantly different veloc- 
ities, where the elongation of the X-ray emission is rotated from 
the axis connecting the two dominants and/or with X-ray peaks 
not coinciding with the location of the dominants (e.g. A2255 
by Burns et al. IT9941 RXJ1314-25 by Valtchanov et al. [20021 
A2744 by Boschin et al. I2006I I. On smaller scales, A773 also 
shows minor but intriguing features. The E-group is not really 
aligned with the direction of the elongation of the X-ray emis- 
sion, but rather shows a PA ~ 90 degrees with respect to the 
cluster centre. Moreover, although the main DM peak roughly 
coincides with that indicated by the two dominant galaxies, it 
indeed is offset by about 1'. 

Further support in favour of an advanced phase of merging 
comes from the luminosity segregations. Although both lumi- 
nous and faint galaxies show two peaks in their velocity dis- 
tribution, only the luminous ones show the E-peak in the 2D 
distribution. Indeed, it is not unusual that galaxies of different 
luminosity trace the dynamics of cluster mergers in a different 
way. The first evidence was given by Biviano et al. ( 19961 see 
also Mercurio et al. 2003 ), who found that the two central dom- 
inant galaxies of the Coma cluster are surrounded by luminous 
galaxies accompanied by the two main X-ray peaks, while the 
distribution of faint galaxies tend to form a structure not cen- 
tred on one of the two dominant galaxies but rather coincident 
with a secondary peak detected in X-rays. Therefore, following 
Biviano et al. ( 1996 ), we might speculate that the merging is in 
an advanced phase, where faint galaxies trace the forming struc- 
ture of the cluster, while more luminous galaxies still trace the 
remnants of pre-merging clumps, which could be so dense to sur- 
vive for a long time after the merging (as suggested by numerical 
simulations; see Gonzalez-Casado et al. 119941 ). 

As discussed in the above sections, the subsystems involved 
in the merger are one cluster and one or two groups with a large 
mass ratio from 4:1 to 10:1, depending on which o" v estimate 
we use for the cluster, with a velocity separation that ranges 
2100 - 2500 km s in the cluster rest-frame. Cosmological sim- 
ulations suggest that these high impact velocities might be oc- 
curring in phenomena of cluster merging (e.g. Crone & Geller 
[19951 ). The value of a merging velocity of ~ 3000 km s is 
also predicted by A^-body numerical simulations of Pinkney et 
al. (Q996) and is probably associated with the phase of the core 
crossing in a cluster-cluster merging process. In the case of 
A773, the extremely large LOS component of the velocity in- 
dicates that the LOS is close to the merging axis. 

As for the direction of the merger, the elongations noted in 
the X-ray emission, the galaxy and DM distributions indicate 
that the merger is not occurring entirely along the LOS, but must 
have at least a small transverse component in the plane of the sky. 
However, since the variety of features do not indicate a unique 
direction we conclude that the geometry of the merger is very 
complex. The likely scenario is that we are looking at a multi- 
ple merger where the N-S direction might indicate the direction 
of an older merger event concerning the cluster and the group 



associated with the D2 dominant galaxy. Dl and D2 might be 
the tracers of the two cores destined to oscillate around the mean 
velocity for a long time (i.e. from two to several 10 9 yr as sug- 
gested by numerical simulations; see, for example, Nakamura et 
al. 119951 Faltenbacher et al. 2006). In fact, a high velocity dif- 
ference between the two dominant galaxies is often suggestive 
of an energetic cluster merger (e.g. Burns et al. 1995) and this 
process is thought to be the cause of the formation of dumbbell 
galaxies in a few merging clusters (e.g. Beers et al. |1992t Flores 
et al. 120001 ). The NE by E-SW by W elongation of the X-ray 
emission might indicate the direction of a more recent merger 
event concerning the eastern group (i.e. that centred on the E- 
peak). We are probably seeing this merging after the phase of 
the core passage. This might explain the offset of the centre of 
the X-ray emission, the offset of the DM peak and the small de- 
viation of the eastern group from the initial NE by E-SW by W 
direction. In this case, the eastern group would have shed its gas 
as a result of ram pressure at the entry into the main cluster at 
its SW by W side and, in fact, we do not see any eastern peak 
in X-ray emission (a scenario already suggested by Govoni et al. 
2004). Further support for this scenario comes from the fact that 
emission line galaxies are preferentially located in the SW clus- 
ter region, the starbursts possibly being induced by the infalling 
subcluster (e.g. Bekki |1999l ). 

One expects that cluster mergers will drive shock waves into 
the intracluster gas of the two subclusters (see Sarazin 20021 for 
a review). The Mach number of the shock is Ai-v s /c s , where 
v s is the velocity of the shock and c s is the sound speed in the 
preshock gas. In the stationary regime we can assume that v s is 
the merger velocity, i.e. v s > 2300 km s taking into account 
the projection factor. We obtain c s < 1 134 km s from the ther- 
mal velocity, i.e. from the observed Tx converted into velocity 
using the usual equipartition of energy density between gas and 
galaxies (see Sect. 5) and taking into account that the observed 
Tx might be enhanced with respect to pre-merging values. We 
therefore estimate M > 2.0. According to our merging scenario, 
this value of M would make front or bow shocks possibly de- 
tectable in X-ray observations. However, in the case of the LOS 
being very close to the merger axis, a front shock would expand 
in an almost parallel direction to the LOS. So because of the 
geometry of the problem we can only observe a global enhance- 
ment of the X-ray temperature in the whole cluster. To date, the 
observational picture is inconclusive. Govoni et al. (2004 ) found 
a very hot region in the western cluster region, but, analysing a 
slightly larger and more homogeneous amount of data, we find 
no significant variation in the X-ray temperature. 

As for the comparison with the radio halo feature, its emis- 
sion is somewhat displaced towards the E-peak (see Fig. [2] and 
Govoni et al. 2004). This is a suggestive indication in favour of 
the connection between cluster mergers and radio halos. Recent 
calculations (Cassano & Brunetti 120051 ) show that a fraction 
~10% of the thermal-cluster energy may be channeled in the 
form of turbulence in case of main (i.e. mass ratio< 3:1) cluster 
mergers and this may power up radio halos. However, a typical 
massive cluster is formed by a series of minor mergers and a 
minor merger may still power the particle acceleration process 
especially if they happen in clusters which are already dynam- 
ically disturbed by precedent mergers (see Fig. 2 and the time- 
evolution of the acceleration coefficient in Fig. 5 in Cassano & 
Brunetti 2005 ). This is probably the case of A773. 
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5.3. LSS structure 

We also analyse the large scale structure (LSS) surrounding 
A773 since the merging direction was shown to be frequently 
correlated with the supercluster environment (e.g. Durret et al. 
1998; Arnaud et al. 120001 1. We find three Abell clusters within 
a circle of 2 degrees centred on A773: A782 of richness class 
R = 2 and A746 and A793 of richness class R — 1 . The clos- 
est one, A782 at 40.6', has five galaxies with known redshift at 
0.20-021 within a circle of 5' as recovered by using we 
therefore assume a cluster redshift of 0.21. Assuming the respec- 
tive redshifts, the distance of A782 to A773 is ~ 40 Mpc thus 
comparable to the typical sizes of superclusters (e.g. Batuski & 
Burns [19831 Zucca et al. [19931 . The direction of A782 is indi- 
cated in Figure [2] the result is that the direction of the elonga- 
tion of X-ray emission and the direction indicated by the eastern 
group are ~ 10-20 degrees different with the axis A773/A782. 
Therefore, we speculate that A773 and A782 lie along the same 
LSS filament where the eastern group also originated. In this 
case, to take into account the geometry of the problem and the 
relative velocity difference of the eastern group and the main 
system, the eastern group lies in front of A773 and has already 
reached the turnaround point, with a second centre passage no 
incipient. 

6. Summary and conclusions 

We present the results of a dynamical analysis of the rich, X-ray 
luminous and hot cluster of galaxies A773 containing a diffuse 
radio halo. 

Our analysis is based on new redshift data for 142 galaxies, 
measured from 107 spectra obtained at the TNG and 37 spectra 
recovered from the CFHT archive in a cluster region within a 
radius of ~ 1.8 Mpc from the cluster centre. We also use 
new photometric data obtained at the INT in a field larger than 
30'x30' . 

We select 100 cluster members with 62250 < v < 68850 
km s _1 and compute a global LOS velocity dispersion of galax- 
ies, cr v = 1394+^ km s _1 . 

The 2D distribution shows two significant peaks separated 
by ~2' along the EW direction with the main, western one lo- 
cated close to the position of the two dominant galaxies and the 
X-ray peak. The velocity distribution of cluster galaxies shows 
two peaks at v ~65000 and ~67500 km s , corresponding to the 
velocities of the two dominant galaxies. The low velocity struc- 
ture has high velocity dispersion - cr v = 800 - 1100 km s _1 - 
and its galaxies are centred on the western 2D peak. The high ve- 
locity structure has intermediate velocity dispersion - <x v ~ 500 
km s _1 - and is characterized by a complex 2D structure with 
a component centred on the western 2D peak and a component 
centred on the eastern 2D peak, these components probably be- 
ing two independent groups. 

Other evidence that A773 is dynamically disturbed comes 
from 1) the presence of two X-ray peaks in the X-ray emission 
in the central cluster region, offset with respect to possible op- 
tical counterparts; 2) the elongation of X-ray emission (PA~80 
degrees) along the NE by E-SW by W direction, while the two 
dominant galaxies indicate a N-S direction; 3) the preferential 
SW location of emission line galaxies. 

According to our likely scenario the high velocity group 
surrounding the high velocity dominant galaxies traces an old 

5 The NASA/IPAC Extragalactic Database (NED) is operated by the 
Jet Propulsion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Administration. 



merger of the main cluster with a group along the N-S direction, 
while the high velocity group surrounding the eastern 2D peak 
of galaxies indicates a more recent merger, but already enough 
advanced, along the NE by E-SW by W direction. Acting on an 
already dynamically disturbed cluster, this second merger event 
would be the likely cause of a very complex observational pic- 
ture and of the radio halo, which is somewhat displaced towards 
the eastern group. 

For the whole galaxy complex we estimate a virial mass 
M = 1.2-2.7 xlO 15 h~l M Q and M(< 1 Mpc)=(5.9- 
ll.l)xl0 14 h~l M Q , values comparable to those of other rich 
clusters at nearby or moderate redshifts. 
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Table 3. Velocity catalog of 142 spectroscopically measured 
galaxies. In Column 1, IDs in italics indicate non-cluster galax- 
ies. IDs 59 and 60 indicate the two dominant galaxies (Dl and 
D2, respectively). 



Table 3. Continued. 



ID a, (5(J2000) 


B 


R 


v Av SC 


09*, +51° 






(km s" 1 ) 



/ 


17 00.60 


42 52.6 


19.77 


18.36 


34762 


70 


e(c) 


2 


17 06.58 


44 07.1 


20.54 


18.86 


64233 


105 




3 


17 13.85 


43 17.0 


22.05 


19.88 


132256 


93 




4 


17 14.04 


43 38.6 


21.24 


19.29 


63970 


99 


k + a 


5 


17 14.50 


43 28.2 


21.73 


20.46 


65568 


123 




6 


17 17.21 


44 43.4 


22.72 


19.95 


132667 


76 




7 


17 19.80 


38 48.8 


21.38 


19.39 


68615 


147 




8 


17 20.86 


38 45.2 


19.63 


18.19 


68559 


94 


e(a) 


9 


17 22.99 


41 18.6 


22.50 


20.17 


64925 


236 




10 


17 23.14 


40 51.6 


22.46 


20.09 


65517 


144 




11 


17 24.05 


38 15.4 


20.67 


19.29 


64503 


87 


e(a?) 


12 


17 26.59 


43 40.4 


21.00 


19.47 


62362 


124 




13 


17 27.86 


37 33.2 


20.53 


18.27 


94246 


76 




14 


17 28.73 


42 23.0 


21.30 


19.28 


62282 


58 


k 


15 


17 28.75 


43 50.2 


21.39 


20.33 


63859 


86 


e(b) 


16 


17 29.71 


42 59.4 


20.93 


18.76 


94885 


76 


e(c) 


17 


17 29.78 


43 08.0 


21.27 


19.32 


68655 


118 




18 


17 31.66 


44 43.1 


21.81 


20.34 


68206 


74 




19 


17 32.50 


42 19.9 


21.77 


19.25 


65171 


111 


e(c?) 


20 


17 32.57 


42 22.2 


22.21 


20.72 


64575 


136 




21 


17 34.25 


41 13.2 


20.64 


18.50 


65848 


69 


k 


22 


17 35.02 


44 32.6 


22.06 


19.73 


67320 


68 




23 


17 35.09 


41 18.6 


20.97 


19.13 


65309 


86 


k 


24 


17 35.38 


38 39.8 


20.44 


18.73 


63729 


105 


k + a 


25 


17 35.59 


41 17.9 


20.66 


18.98 


63067 


99 




26 


17 36.10 


44 18.2 


20.75 


19.44 


58492 


81 




27 


17 36.29 


45 27.0 


21.89 


19.73 


67997 


87 




28 


17 37.08 


42 53.6 


22.26 


19.87 


93910 


73 




29 


17 37.37 


42 36.4 


22.85 


20.73 


64940 


212 




30 


17 37.99 


44 17.5 


20.22 


18.08 


66016 


51 


k 


31 


17 38.66 


43 25.7 


20.60 


19.30 


64889 


85 


e(a?) 


32 


17 38.74 


46 58.4 


19.67 


18.84 


20854 


114 


e(a?) 


33 


17 39.22 


46 53.8 


20.93 


18.80 


64216 


118 


k 


34 


17 39.48 


41 49.6 


21.97 


20.20 


67577 


70 




35 


17 40.73 


38 18.6 


19.95 


18.72 


67396 


86 




36 


17 40.80 


41 31.6 


22.54 


20.13 


65289 


206 




37 


17 41.06 


43 58.4 


22.19 


20.68 


124666 


24 


e(c?) 


38 


17 42.79 


43 41.2 


20.90 


19.77 


55811 


93 


e(c) 


39 


17 44.04 


40 21.0 


22.51 


20.52 


67468 


148 


e(c?) 


40 


17 44.66 


43 26.4 


21.73 


19.63 


65758 


70 


k 


41 


17 45.72 


42 55.1 


21.66 


19.52 


64248 


72 


k 


42 


17 46.32 


43 13.1 


21.26 


18.90 


65547 


50 


k 


43 


17 46.85 


45 29.5 


20.74 


18.64 


65902 


54 


k 


44 


17 47.14 


40 59.2 


21.26 


19.03 


65661 


106 


k 


45 


17 48.02 


46 35.4 


20.30 


18.65 


66653 


112 


k 


46 


17 48.22 


46 13.1 


20.85 


18.76 


66848 


60 


k 


47 


17 48.62 


43 43.7 


22.07 


20.81 


64023 


90 




48 


17 48.82 


40 59.9 


19.90 


18.78 


44429 


81 


e(a) 


49 


17 48.89 


44 53.5 


20.05 


18.73 


55538 


92 


e(a) 


50 


17 49.46 


43 08.4 


21.69 


19.18 


68830 


76 


k 



ID 


a,<5(J2000) 


B 


R 


V 


Av 


SC 




09*, +51° 






(km 






51 


17 AQ f,f, 


zl1 Ozl 

H 1 Ut.7 


19.88 


17.64 


65088 


62 


k 


52 


1 7 4Q 89 

1 / H77.0Z 


zn 5i n 

HJ Jl.u 


22.67 


21.27 


67857 


58 e(c?) 


53 


1 7 5f) 54 


zn 96 zl 

HJ ZU.4 


20.95 


18.63 


63900 


38 


k 


54 


1751 fin 


zn 1 6 n 

HJ 1U.U 


22.27 


19.68 


63696 


64 




55 


17 51 Q6 

1 / J 1.7U 


zl1 XX 7 

HI JO. / 


20.13 


19.04 


114971 


75 


e(a) 


56 


17 59 56 


zn 50 ? 

HO JU.Z 


20.84 


18.31 


64849 


46 


k 


57 


1 7 5^ OQ 

1 / JJ.U7 


ZL9 55 ZL 
HZ JJ.H 


20.33 


18.05 


67364 


72 


k 


58 


17 5^ 16 

1 / JJ.1U 


44 08 

HH U0.7 


20.81 


18.59 


64194 


74 


k 


59 


1 7 ca 9ft 


zn "Xf. 5 
Hj jD.J 


18.82 


16.49 


65131 


50 


k 


60 


17 5^ ^ 


zlzt Of) 6 

HH UU.U 


18.93 


16.73 


67765 


56 


k 


61 


17 5^! ^8 


zLzl 94 

HH ZH.U 


21.40 


19.05 


65152 


124 


k 


62 


17 54 17 

1 / JH. 1 / 


zn 97 5 


21.26 


18.87 


67578 


126 


k 


63 


17 54 17 


zn 55 6 


19.90 


17.83 


66477 


28 


k 


64 


17 54.24 


44 52.8 


20.44 


18.58 


68485 


66 


k 


65 


1 7 54 89 


zL9 zl6 8 

HZ HU.O 


22.46 


19.99 


79178 


158 




66 


1 7 54 84 


4.7 0^ 8 

H/ UJ.O 


18.84 


17.77 


21910 


116 


e(a) 


67 


17 55 54 


zn 5zi i 

H J Jt. 1 


20.77 


18.49 


66543 


30 


k 


68 


1 7 5 5 Q7 

1 / JJ.7 / 


4.1 51 

HI Jl.U 


20.85 


19.24 


55993 


87 


e(c) 


69 


17 56 76 
1 / ju. / u 


zn 98 6 
tj zo.u 


21.13 


18.89 


64437 


57 


k 


70 


17 57 1f) 

i / j i . iu 


4.6 1f) Q 

HU 1U.7 


21.20 


19.00 


68250 


50 


k 


71 


17 57 89 
1 / J 1 .oz 


zn 97 zl 

HO z I .H 


23.41 


21.39 


66885 


130 




72 


17 58 9fi 


4.1 ^5 5 

HI JJ.J 


21.01 


18.59 


67420 


39 


k 


73 


1 7 5Q 76 

1 / J7. / U 


4.0 0^ 7 
HU vtj. 1 


21.76 


19.54 


67738 


64 


k 


74 


i c on 1 zi 

1 UU. 14 


46 TX 9 

HU ZJ.Z 


21.24 


19.12 


68085 


51 


k 


75 


1 o fin 94 

1 UU.ZH 


44. 1 1 n 

HH ll.U 


20.89 


19.02 


62377 


46 


k + a 


76 


i s on 60 
1 o uu.ou 


44 07 4 
HH U / .H 


21.30 


18.86 


64978 


33 


k 


77 


18 01 97 
10 U 1 .z / 


44 1 ^ 6 

HH lj . U 


21.58 


19.15 


68682 


69 


k 


78 


18 01 ^7 

1 O U 1 . J I 


49 95 Q 

HZ ZJ.7 


21.54 


19.24 


66377 


69 


k 


79 


18 01 75 


40 1.A 

HU JH.U 


21.27 


18.85 


64829 


72 


k 


80 


1 8 0? 86 
1 uz.ou 


40 07 1 

HU U / . J 


21.58 


19.55 


63748 


87 


k 


81 


1 8 0? 86 
10 uz.ou 


41 4f) 6 

HI HU . U 


20.27 


18.78 


68358 


44 


a + k 


82 


1 8 O^ ZL6 
10 Uj.tD 


4^ 51 9 
Hj Jl.z 


20.20 


17.98 


67742 


23 


k 


83 


1 8 rn 55 

1 O UJ.JJ 


4^ ^n 9 

tj JJ.Z 


19.95 


17.77 


67969 


54 


k 


84 


1 8 O^ 77 
1 Uj. / / 


4^ 98 9 
Hj ZO.Z 


21.29 


19.24 


66659 


166 




85 


1 8 Ozt 1 

10 Irt.Ul 


49 54 f) 

HZ JH.U 


21.50 


19.41 


65774 


66 


k 


86 


1 8 OA 1 O 
1 UH. 1U 


4^ 1 9 7 
Hj 1Z. / 


21.05 


18.76 


64671 


46 


k 


87 


1 8 Ozt zl? 

lO UH.HZ 


AT. AA A 
tj HH.H 


19.59 


17.53 


67673 


36 


k 


88 


1 8 05 06 
1 UJ.UO 


AT, AQ 8 

4J Hy.o 


20.04 


17.98 


67047 


24 


k 


89 


18 05.38 


41 04.2 


21.80 


20.04 


64928 


112 


e(a) 


90 


18 05.57 


46 53.4 


20.50 


18.51 


67098 


52 


k + a 


91 


18 05.64 


39 59.0 


21.46 


19.38 


12022 


72 




92 


18 05.71 


45 18.0 


20.50 


18.52 


66287 


66 


e(c) 


93 


18 05.78 


43 40.1 


20.22 


18.20 


65414 


24 


k 


94 


18 06.22 


42 00.0 


20.90 


18.91 


66502 


60 


k 


95 


18 07.08 


44 17.2 


20.51 


18.23 


65081 


27 


k 


96 


18 07.27 


49 13.8 


19.56 


18.44 


44491 


88 


e(b) 


97 


18 08.57 


43 42.6 


21.16 


18.98 


67860 


28 


k 


98 


18 08.57 


48 51.5 


20.36 


19.28 


44498 


70 


e(c) 


99 


18 08.64 


48 11.3 


0.0 


0.0 


65850 


70 


k 


100 


18 09.79 


41 48.1 


20.36 


18.70 


67183 


69 


e(a) 
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Table 3. Continued. 



ID a,<5(J2000) 


B 


R 


v Av SC 


09*, +51° 






(km s" 1 ) 



101 18 

102 18 

103 18 

104 18 

105 18 

106 18 

107 18 

108 18 

109 18 

110 18 

111 18 

112 18 

113 18 

114 18 

115 18 

116 18 

117 18 

118 18 

119 18 

120 18 

121 18 

122 18 

123 18 

124 18 

125 18 

126 18 

127 18 

128 18 

129 18 

130 18 

131 18 

132 18 

133 18 

134 18 

135 18 

136 18 

137 18 

138 18 

139 18 

140 18 

141 18 

142 18 



10.13 
10.54 
10.61 
11.45 
12.37 
12.98 
14.38 
14.54 
14.81 
15.17 
15.84 
15.98 
17.04 
17.11 
18.34 
18.55 
19.66 
21.00 
22.13 
22.63 
23.18 
23.42 
23.47 
23.74 
23.77 
23.93 
24.57 
24.72 
25.10 
27.72 
27.72 
27.73 
29.57 
30.21 
30.46 
31.73 
35.30 
40.66 
41.33 
44.45 
44.86 
49.37 



,45 18.0 
, 42 26.3 
,43 14.5 
,43 44.8 
,47 48.1 
, 42 49.0 
,46 43.7 
,43 20.9 
,43 49.4 
,43 21.8 
,40 44.4 
,45 33.5 
, 42 50.0 
,43 28.2 
,47 53.5 
, 43 49.0 
, 42 48.2 
,43 05.5 
,47 56.0 
,43 09.5 
,48 33.1 
,45 22.0 
, 45 56.9 
,43 30.4 
, 47 24.9 
,47 19.0 
,47 54.5 
,42 01.4 
,49 09.1 
,41 13.2 
,46 05.9 
,44 27.7 
, 40 42.6 
,47 21.9 
,44 23.6 
, 42 49.3 
',45 26.6 
,44 12.1 
,43 10.6 
,44 14.3 
,43 09.8 
,45 19.8 



21.29 
20.80 
21.33 
20.47 
21.37 
21.57 
20.83 
20.27 
21.38 
19.77 
20.32 
20.77 
20.07 
22.02 
21.75 
20.23 
22.30 
20.67 
21.81 
21.79 
20.79 
21.05 
21.47 
22.73 
20.20 
21.15 
21.50 
21.84 
20.78 
22.15 
20.52 
21.36 
21.57 
22.38 
20.86 
20.60 
19.77 
20.34 
21.11 
21.78 
19.92 
20.72 



19.12 
18.79 
19.41 
19.60 
19.81 
19.90 
19.04 
18.10 
19.22 
17.74 
18.87 
19.23 
18.04 
19.67 
19.39 
17.59 
19.84 
18.38 
19.84 
20.76 
18.47 
18.84 
19.25 
20.62 
19.38 
19.34 
20.20 
20.19 
19.43 
20.43 
19.52 
19.38 
19.48 
20.21 
19.64 
18.51 
17.60 
18.34 
20.13 
19.71 
18.21 
19.30 



65867 

66845 

67062 

55929 

101181 

63437 

63833 

55650 

64637 

55648 

55650 

47297 

55700 

78329 

64910 

86619 

65556 

67354 

64686 

44170 

64095 

64388 

74424 

67352 

58662 

63302 

65857 

102297 

128313 

130562 

21763 

65252 

67741 

138204 

52800 

64282 

78472 

52942 

29057 

65279 

52748 

113177 



147 

24 
69 
80 
88 
96 
84 
34 
69 
60 
84 
102 
38 
64 
60 
25 
36 
35 
117 
159 
57 
69 
60 
50 
141 
86 
92 
47 
96 
72 
130 
110 
33 
74 
90 
66 
56 
64 
104 
122 
63 
231 



k 

k 
e(b) 
e(c?) 

k + a 

k 

k 
e(a) 

k 



k 
k 
k 
k 
k 
e(b) 
k 
k 
k 

e(b) 
k 
k 
k 

e(c?) 



k 
k 

e(c) 
k 

k + a 
k 

e(c) 

k 
e(a) 



